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Abstract. At the distance of NGC 5128 (3.6±0.2 Mpc) 
it is possible to resolve globular clusters with high resolu- 
tion imaging from the ground, thus allowing the globular 
cluster candidate selection primarily through their mor- 
phological properties. I report the discovery of 71 globu- 
lar clusters in NGC 5128 on VLT UTl+FORSl images, 
including the faintest members (My —5) known to 
date in this galaxy as well as 5 previously known clus- 
ters. U- and V-band photometry has been measured for 
all the candidates and the luminosity function, spanning 
-10.1 < Mv < -4.9 and -9.3 < Mu < -3.3, con- 
structed. These are the deepest globular cluster luminos- 
ity functions in an elliptical galaxy determined so far. 
The Kolmogorov-Smirnov statistics show that the differ- 
ence between the globular cluster luminosity functions of 
NGC 5128 and the MW is not larger than the difference 
between the ones of M31 and the MW. The (U-V)o color 
histogram shows a bimodal distribution. For 23 globular 
clusters I obtained K-band images with SOFI at the NTT 
in La Silla. Their positions in the (U— V) vs. (V— K) color- 
color diagram indicate that they are indeed old globular 
clusters. Assuming that the globular clusters in NGC 5128 
span a similar age range as the ones in the Milky Way and 
adopting a linear fit between the metallicity and (U— V)o 
color, the metal-rich clusters peak at [Fe/H]= —0.6 dex 
and the metal-poor ones peak at [Fe/H]= —1.7 dex. 
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1. Introduction 

Globular Clusters are among the oldest objects in the Uni- 
verse. Moreover, they are intrinsically luminous and found 
around almost all galaxies in which searches have been 
made. There are growing evidences for globular cluster for- 
mation during major star formation episodes (Schweizer 



1997; Larsen & Richtler 1999). This makes them useful 



tracers of the formation and the evolution of galaxies (Har- 



ris |1991| , pOOOj ; Kissler-Patig et al. |1998t van den Bergh 
^OOOl) . 



The globular cluster luminosity function (GCLF) of 
a galaxy is the relative number of its globular clusters 
per luminosity (or magnitude) interval. It is usually uni- 
modal and nearly symmetric with a peak at a characteris- 
tic (turn-over) magnitude. Empirically it was found that 
turn-over magnitudes for galaxies at the same distance 
(e.g. in the same galaxy cluster) differ by < 0.1, which 



makes the G CLF a good standard candle (Harris 1991 
Jacoby et al. 1992| ). However, the full calibration of the 



GCLF is incomplete. The essential problem is that the 
galaxies at large distances, for which the GCLF method is 
mostly used, are giant ellipticals in which globular clusters 
are found in largest numbers and in which there are usu- 
ally fewer problems with extinction, whilst the primary 
calibration of the standard candle My rests on the MW 
and M31, which are spiral galaxies. In fact, evidence that 
GCLFs depend on Hubble type and environment has been 



claimed (Fleming et al. 1995), but the former is proba- 
bly due to differences in mean metallicity of the globular 
clusters and the latter was never reliably demonstrated 



(Kissler-Patig |2000| ). 

There is strong observational evidence that most lu- 
minous giant ellipticals have a bimodal globular cluster 
metallicity distribution (Gebhardt & Kissler-Patig 1999; 
Kundu 1999). The globular cluster specific frequencies 



(the number of clusters norm alized to the galaxy lumi- 
nosity; Harris & van den Bergh [l981 ) ar e typ ically higher 
in ellipticals than in spirals (Harris |l 991 2000). To explain 
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the high number of globular clusters and the presence of 
several populations indicated by the bimodal metallicity 
distributions, different globular cluster formation mech- 



anisms have been proposed (e.g. Ashman & Zepf 1992; 
Cote et al. 1998| ; Forbes et al. 1997). The number, chemi- 
cal abundance and spatial distribution of a globular clus- 
ter system (GCS) put strong constraints on the formation 
of its parent galaxy. 



At th e distance of 3.6 Mpc (Soria et al. 199f ; Harris et 
al. 1999) NGC 5128 is the closest giant elliptical galaxy. 
The first globular cluster in NGC 5128 was discovered by 
Graham & Phillips (1980). Its brightness was high with 
V = 17.2 ± 0.1 mag, which corresponds to My = —10.6 
mag. Subsequent investigation of photographic plates (van 
den Bergh et al. |l98l| , Hesser et al. |1984| , |1986D revealed 
a rich GCS in this galaxy. The estimated number of glob- 
ular clusters in NGC 5128 is N 
compared with N- 



-1700 (Harris et 
-150 for the MW. 



1984) 



On the basis of the near IR colors for a magnitude- 
limited sample of 12 globular clusters, Frogel (198C, 1984) 
found one third of the sample to be super metal rich 
(above solar metallicity) and luminous. However, the 
Washington photometry of 62 spectroscopically confirmed 
clusters (Harris et al. 1992) as well as a spectrophotom- 
etry for 5 clusters (Jablonka et al. 1996[ ) showed less ex- 
treme metallicity. Harris et al. (1992) found a large spread 
in metallicities with a mean of [Fe/Hjc-Ti = —0.8 ± 0.2 
dex. They also find evidence for a weak spatial metallicity 
gradient, but state that the same could be produced by 
unknown reddening effects. Zepf & Ashman ( 1993| ), using 
the same photometric data, suggest that the metallicity 
distribution of globular clusters in NGC 5128 is bimodal, 
with the higher metallicity peak at [Fe/H]=-|-0.25 dex, and 
proposed that this was due to the formation of a second 
population of globular clusters in a previous merger which 
created NGC 5128. 

Several photometric studies in the inner regions of 
NGC 5128 have been conducted, with the aim of con- 
straining the dependence of globular cluster metallicity on 
galactocentric radius and of finding super-metal-rich clus- 
ters (Minniti et al. [1996| ; Holland et al. \L99^ . These stud- 
ies suggest the existence of radial gradients and bimodal 
distribution in globular cluster metallicity. However, their 
conclusions are limited due to the restricted samples, pos- 
sible reddening within NGC 5128 that is especially severe 
in the central parts around the prominent dust lane, and 
possible contamination by stars within NGC 5128 as well 
as background galaxies. 

In this paper, I present a sample of globular clusters 
from two fields in the halo of NGC 5128. Deep images 
combined with high resolving power of the VLT and the 
proximity of NGC 5128 allow the identification of globu- 
lar clusters primarily on the basis of their non-stellar PSF. 
Constructing the GCLF for a statistically significant num- 
ber of clusters that span the whole range of magnitudes. 



I pursue the question of the uniqueness of GCLF in ellip- 
tical galaxies. 

This paper is organized as follows. The data, the pho- 
tometric calibrations and completeness tests are described 
in Sect. ||. In Sect. ^ the color-magnitude diagram for the 
NGC 5128 globular clusters is compared to the ones of 
the MW and M31. The luminosity function is presented in 
Sect. ^ and the color distribution in Sect. |[ In the latter, 
the mean abundances are determined for the metal-poor 
and metal-rich populations. The results are summarized 
in Sect. ^. 

2. The Data 

2.1. Observations and Data Reduction 

I used the VLT Antu (UT1)+F0RS1 (FOcal Reducer/low 
dispersion Spectrograph) to obtain the Bessel U- and V- 
band photometry of globular clusters and stars in two dif- 
ferent fields in NGC 5128. The field of view was 6'.8 x 6.'8 
and the pixel scale O'.'2/pixel. The FORSl detector is a 
2048x2048 CCD with 24^m pixels. 

Field 1 was centered on the prominent N-E shell ~14' 
away from the center of the galaxy, while Field 2 was cho- 
sen to overlap with Soria et al. ( 1996 ) HST observations 
at a distance of ^^9' from the center of the galaxy. Ob- 
servations were carried out on the 11th and 12th of July 
1999 in service mode. Pairs of images of 15 min exposure 
time were taken per field in each filter. The journal of 
observations is given in Tab. |. 

In addition to VLT data, I obtained a 45 min long ex- 
posure of Field 2 in using SOFI at the New Technology 
Telescope (NTT) at the ESO/La Silla Observatory. The 
field of view of SOFI is 4^94 x 4f94 and the pixel scale 
O'.'292/pixel. This field is smaller than the one of FORSl 
and so only about 40% of objects in one field have Ks-band 
magnitudes. 

Both fields were also observed with the Wide Field Im- 
ager (WFI) at the 2.2m telescope on La Silla. The data 
from WFI observations are used here to assess the V-band 
magnitude of some of the brightest cluster candidates that 
were saturated on V-band FORSl images. I calibrated 
the WFI photometry by comparing magnitudes of globu- 
lar clusters well exposed (but not saturated) on VLT and 
WFI images in order to determine the zero point of WFI 
observations. 

For optical data (VLT and 2.2m images), the standard 
image processing, including the overscan and bias subtrac- 
tion and the flat-field correction, was performed within the 
IRAF[| environment. The pairs of exposures in each filter 
were registered using the imalign task and averaged to 
obtain the final images. Cosmic rays were rejected in this 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under contract with the 
National Science Foundation 
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Table 1. Journal of Observations 
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(h min sec) 
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dd/mm/yy 
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(sec) 
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process using the crreject algorithm within the imcombine 
task in IRAF. Bad pixels were masked out. 

The standard procedure in reducing infrared (IR) data 
consists of sky subtraction, flat-field correction, register- 
ing and combining the images. Good sky subtraction in 
a crowded field like that of a galactic halo is particularly 
important. For that step I used the DIMSUM package 



(Stanford et al. |1995| ) within IRAF. In DIMSUM the sky 
subtraction is made in two passes. In the first one a median 
sky is computed for each image from 6 temporarily closest 
frames. The shifts between the sky subtracted frames are 
then computed and all the images stacked together using 
a rejection algorithm to remove cosmic rays. An object 
mask is computed for the coadded image and then shifted 
back in order to create object masks for the individual 
frames. In the second pass, the sky subtraction is made 
using the object masks to avoid overestimation of the sky 
level. These masks are also used to check that the bright 
object cores were not removed as cosmic rays in the previ- 
ous pass. After the mask-pass sky subtraction, the frames 
are registered with imalign and combined with imcombine 
task in IRAF. 

2.2. Cluster identification and photometry 

NGC 5128 is close enough that its globular clusters can be 
distinguished from stars on the basis of their slightly non- 
stellar appearance on high resolution ground-based images 
taken in excellent seeing conditions. At the distance of 
NGC 5128 (3.6 Mpc) a typical MW globular cluster with 
a mean King core radius of Tc = 2.3 pc had, on the best 
seeing VLT images, a FWHM of 2.6 pix, which is ~ 0.4 
pix bigger than the stellar FWHM. The excellent seeing, 
especially on the V-band images of Field 2, allowed me 
to identify globular cluster candidates with Tc as small 
as 1.7 pc, corresponding to observed FWHM of 2.4 pix 
(see f2.GC-6 in Tab. purely on the basis of their non- 
stellar appearance through the following procedure within 
IRAF: (i) first the daofind algorithm was run in order to 
detect all the stars and clusters in the image; (ii) pho- 
tometry through a 3 pixels aperture was performed on 
all the objects; (iii) relatively bright, isolated, stellar ob- 
jects were used to create the PSF, which was then fitted 
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Fig. 1. Residual image of a 51" x 51" portion of Field 2. 
On the left panels the V-band image (top) and its residual 
image (bottom) and on the right panels the corresponding 
U-band images are shown. Three globular cluster candi- 
dates are circled. The left-most globular cluster is slightly 
overexposed in the V-band image. Its classification as a 
globular cluster follows from the U-band residuals. 



to all the objects, subtracted and the residual images cre- 
ated. DAOPHOT II (Stetson |l987| ) permits creation of the 
variable PSF, which was necessary for the U-band images 
that showed fixed pattern residuals upon subtraction of 
the non-variable PSF. In order to create a variable PSF I 
used 35 stars per image. 

Because of their slightly larger FWHM and non-stellar 
PSF, globular clusters were easily detected on residual im- 
ages created by subtracting the objects fitted with a stel- 
lar PSF. This subtraction was performed using the allstar 
task. Globular cluster residuals were oversubtracted in the 
center and undersubtracted in the wings (Fig. |l|). Images 
were visually inspected and the residuals in the U- and 
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V-band images compared. This was necessary to discard 
obvious galaxies from the candidate Ust. 

After visual inspection of the residual images, 50 and 
81 objects in Field 1 and 2, respectively, were retained 
as candidates for globular clusters. Good globular cluster 
candidates, the ones that have round ring-like residuals, 
are listed in Tab. |[ Other extended objects for which 
the visual distinction was not very clear, some of them 
being blended stars or round galaxies, but some of which 
might be globular clusters as well, are listed in Ta b. ^. 
The Sextractor programme (Bertin & Arnouts |l996 ) was 
used to measure precise coordinates, FWHM {fu and fv 
columns in Tab. ^ and ^ and ellipticities {ey and ejj 
columns) for these objects. 

I performed aperture photometry within IRAF for all 
the objects. Since the fields were quite crowded, stellar 
objects had to be removed in order to obtain more pre- 
cise photometry of globular clusters. To do so, 1 first sub- 
tracted all the stars (this time excluding the globular clus- 
ter candidates from the aperture photometry list) from 
the original images using the allstar task. Aperture pho- 
tometry of globular cluster candidates was then measured 
through circular apertures of 16, 20 and 25 pixel radii. 
In some cases a bright saturated star or another globular 
cluster candidate was found near the cluster, or the latter 
was located near the edge of the image and thus only a 3 
pix aperture magnitude was used. The 3 pix magnitudes 
were later extrapolated to 25 pix magnitude values using 
an aperture correction calculated as a median correction 
for more isolated globular clusters. The error in magni- 
tude for these clusters is larger and hard to quantify, since 
the aperture correction depends on the concentration and 
compactness of the cluster. The positions and magnitudes 
of globular clusters are presented in Tab. || and Tab. ^. 
The clusters with magnitudes measured through a smaller 
aperture and later corrected are flagged with "b" (near the 
bright star or another globular cluster) or "e" (near the 
edge of the image). 

2.3. Photometric calibration 

For the photometric calibration of the VLT images, a total 
of 14 stars in 4 different Landolt (1992) fields, spanning 
the color range —1.321 < {U — V) < 4.162, were observed 
during the two nights. The following transformations were 
calculated: 

Uinst = ?7- 24.262(±0.089)H-0.379(±0.068)*X 

-0.042(±0.007) * (?/ - y) (1) 



V ~ 27.348(±0.042) + 0.213(±0.034) * X 



(2) 



where X is the mean airmass of the observations, Umst and 
Vinst are the instrumental magnitudes and U and V mag- 
nitudes from the Landolt ( 1992| ) catalogue. The one sigma 
scatter around the mean was 0.031 mag for the U band and 
0.023 mag for the V band (Fig. ^ left panel). Adding the 
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Fig. 2. Photometric calibration of Landolt ( 1992| ) stars 
during the nights of observations in U and V filters. Left 
panels: the scatter for the calibration without the color 
term; right panels: the scatter for the calibration with the 
color term. 



ZPT = 22.34 ±0.01 
k - -0.058 ±0.007 
a= 0.028 
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Fig. 3. Photometric calibration of Persson et al. ( 1998 ) 
stars in K, filter. 



color term (U— V) in the transformations slightly reduces 
the scatter in the V-band to 0.017 mag (Fig. || right panel). 
However, the calibration equation without the color term 
for the V-band was prefered, since the data in that filter 
go much deeper and U magnitudes for some objects could 
not be measured accurately enough. 

SOFl images were calibrated using 6 standards from 
the Persson et al. ( 199^ ) list observed during the night and 
spanning an airmass range from 1.1 to 3.4. Each standard 
star was observed at 5 different positions on the IR-array. 
In this way total of 30 independent measurements was 
obtained. The least square fit (Fig. ^) yielded the following 
calibration equation: 



K, - 22.34(±0.01) + 0.058(±0.007) * X (3) 
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2.4- Completeness 



Table 2. Comparison with the data from literature. 



I made completeness tests for detection of stars on the U- 
and V-band images. Using the IRAF task addstar within 
the DAOPHOT package I added 400 stars each time to 
each of the U- and V-band images, all with the same 
magnitude and uniformly distributed over the whole field, 
and re-computed the photometry. The completeness and 
the accuracy of photometry were measured at 0.5 mag 
bins. The 90% completeness for stellar sources is reached 
at magnitude 25 in the V-band and 23.5 in U-band. The 
completeness drops to 50% at V~ 26 and U~ 24.5. More 
detailed results of completeness tests will be presented in 
the forthcoming paper where the stellar photometry of 
NGC 5128 is analyzed (Rejkuba et al. |2001b| ). 

Globular clusters have larger PSFs than stars, their se- 
lection is based mainly on the recognition of a "globular- 
like" residual on a subtracted image and thus their iden- 
tification as such is more subjective. In order to check the 
completeness level of globular cluster detections, I made 
similar simulations as for the stars. The fainter globular 
clusters, that are of interest in simulation, are not resolved 
and thus I just treat them as stars. To add a correct, 
typical globular cluster PSF, the best cluster candidates' 
PSFs were used to create simulated images. In order to 
investigate the dependence of detection on core radius, 
simulations with different typical PSFs were made. Two 
sets of simulations were done using the PSFs of clusters 
with small core radii, fl.GC-17 and f2.GC-4, to add ob- 
jects, while in the other two sets much less concentrated 
clusters fl.GC-10 and f2.GC-36 were used. At each pass 
100 clusters were added, distributed uniformly all with the 
same magnitude. The completeness of the globular cluster 
detection was measured for each half magnitude. 

The detection depends quite significantly on core ra- 
dius in the sense that more small Vc clusters are lost at 
faint end. The incompleteness magnitude is defined as 
magnitude where more than 50% of the added clusters are 
not detected as such (they are actually all detected, but 
their residuals are not recognized as globular cluster-like, 
because they are lost in noise). It sets in at V'^22.7 mag for 
clusters like fl.GC-17 and f2.GC-4, but most of the larger 
clusters, like fl.GC-10 and f2.GC-36 can still be detected 
at 90% level at magnitude V=23. The incompleteness for 
clusters with large Tc sets in at magnitude V=23.5. In the 
simulations clusters could be detected with magnitudes as 
faint as V=24. However, at those magnitude limits it is 
almost impossible to distinguish globular cluster residuals 
from galaxies. 

From the simulations I conclude that clusters fainter 
than 23 mag would probably not be detected. MW 
clusters corresponding to these magnitudes are Pal 1, 
Pal 13, Terzan 1, E 3 and AM 4. However, at magnitude 
V = 22.0 (~ 2 magnitudes past the GCLF peak) more 
than 90% of the clusters are detected. 



IDl 


ID2 


Bpg 




Vwash 


VwFI 


(U-V) 


fl.GC-13 


C36 


19.04 


18.24 


18.27 


18.39 


0.91 


fl.GC-16 


C37 


19.24 


18.17 


18.46 


18.47 


1.33 


fl.GC-23 


C38 


19.00 






18.40 


1.16 


f2.GC-26 


C15 


19.76 


18.62 




18.58 


1.40 


f2.GC-81 


C12 


18.82 


17.74 


17.79 


18.01 


1.47 



2.5. Comparison with published data 

The most recen t studies of the GCSin NGC 5128 (Hol- 
land et al. 1999, Alonso & Minniti 1997) have investigated 



the central regions of the galaxy away from the fields stud- 
ied here. The only photometry of globular clusters found 
in common with the data presented here is the photo- 
gra phic p l ate m easurements of five clusters in Hesser et 
al. ( 1986 , 1984) and the Washington p hotom etry of the 
three of these 5 clusters in Harris et al. (1992). The clus- 
ters in common are C36 , C37 , C38 , C12 and C15 (num- 
bering from Hesser et al. 1984, 1986 ). While in the Hesser 
et al. papers only the (B) magnitude is given, Harris et 
al. (|1992| ) report also photographic V-band magnitudes as 
well as Washington photometry for C12, C36 and C37. In 
order to compare the Washington photometry with Bessell 
V magnitudes, I used the transformation equation (Harris 
et al. |1992| ): V = Ti + 0.66(M - Ti). 

The clusters in common are among the brightest ones 
and in this study their magnitudes were measured on 
2.2m-|-WFI images. The comparison between the data 
from the literature and my measurements is given in 
Tab. H. In columns 1 and 2 the identification of the objects 
in common is given, Bpg magnitudes are from Hesser et al. 
( 1984 , 1986 ), Vpg and Vwash are photographic magnitudes 
and the V-band brightness derived from Washington pho- 
tometry. The last two columns present my measurements. 
There is a good agreement between the V-band magnitude 
derived from the Washington photometry and my mea- 
surements. The largest difference is for cluster C12 (f2.GC- 
81), which I measured 0.2 magnitudes fainter, probably 
due to the vicinity of a bad column on the WFI chip. 



2.6. Selection of globular clusters 

Typical Galactic globular clusters have mean King core 
radii Tc = 2.3 ± 3.6pc, mean tidal radii rt = 45.1 ± 32.4pc 
and mea n ellip ticities e = 0.07 ± 0.01 (J une 19 99 update 
of Harris |l996| cat alogue ; White & Shawl 119871 ). 

Holland et al. ( p^99| ) used HST WFPC2 data to find 
globular clusters in the central parts of NGC 5128 and 
measured their structural parameters. They did not find 
any significant difference between the MW and NGC 5128 
globular clusters' King core radii (rc), tidal radii (rt) and 
half-mass radii (r^,). They also found no trend in Tc, rt or 
r/i with galactocentric distance, contrary to what one may 
expect due to stronger tidal forces closer to the galactic 
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Table 3. Coordinates, photometry, FWHM (/y and fu), ellipticities (ey and cu) and projected galactocentric distance 
of globular cluster candidates. 
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Table 3. continuation. 
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comment: b = near bright (saturated) star 

e = near edge of the field 

u = marginal or no detection 

V = object saturated in V; V mag from WFI data 

n = object does not satisfy all the criteria for globular cluster classification 
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Fig. 4. FWHM vs. ellipticity: objects retained as good 
globular cluster candidates after the visual inspection of 
the residual images of Field 1 and 2 are plotted with filled 
circles. Empty circles are used for other extended objects. 



center. However, they measured higher mean ellipticities 
for their globular cluster candidates with respect to the 
MW globular clusters. Their samples are restricted only 
to the brighter clusters in NGC 5128, corresponding to 



the absolute magnitude My = —6.4. For comparison, the 
present sample reaches almost 2 magnitudes deeper. 

I measured the FWHM and ellipticity for all the ob- 
jects in Tab. ^ a nd ^ using the Sextractor programme 
(Bertin & Arnouts 1996| ). U-band images had variable PSF 
across the frame and the objects were, in some places, 
elongated, thus showing higher ellipticity than in the V- 
band images. 

The FWHM vs. ellipticity for aU the objects in Fields 1 
and 2 is plotted in Fig. ^. The upper panel shows the 
results obtained from the V-band, the lower panel the ones 
obtained from the U-band. The vertical line indicates the 
mean FWHM for stars. The objects from Tab. |, some 
of which might as well be blended stars or galaxies are 
plotted as open symbols, while the filled dots were used 
for objects from Tab. || (i.e. good candidates). Objects 
that were slightly saturated on the V-band VLT images 
are shown only in the U-band diagram (lower panel). Due 
to the poorer seeing conditions during the observations 
with the 2.2m telescope, their profiles are more similar to 
the stellar profiles on the WFI images. 

In order to classify an object as a globular cluster all 
of the following criteria were imposed: 

1. non-stellar PSF that leaves ring-like residuals on sub- 
tracted images (Fig. |l|); 
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Table 4. Same as Tab. ||, but for other extended objects, some of which might as well be globular clusters. 
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comment: b = near bright (saturated) star 

e = near edge of the field 

u = marginal or no detection 

V = object saturated in V; V mag is from WFI data 
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Fig. 5. Color-Magnitude Diagram for all objects found in 
Field 1 (diamonds) and in Field 2 (triangles). Objects re- 
tained as good globular cluster candidates after the visual 
inspection of images are plotted as filled symbols and the 
rest as open ones. 

2. FWHM>mean stellar FWHM and ellipticity(V)<0.4 
(Fig. I) 

3. < (t/ ~ V)o < 2.5 (Fig. 1) 

The selection of objects with ellipticity measured on the 
V-band images smaller than 0.4 excludes galaxies from the 
sample. The color range in the last criterion corresponds 
to the observed range of the (U— V)o colors for the MW 
and M31 globular clusters. This range does not include the 
reddest M31 clusters (Fig. ||) because they were corrected 
only for the foreground average extinction towards M31 
and thus might still be reddened. The cutoff at (U— V)o — 
2.5 takes into account the bulk of M31 clusters around 
Vq ~ 21. The objects that did not pass all of the above 
criteria were flagged with "n" in Tab. ||. 

3. Color-Magnitude and Color-Color Diagrams 

The UV color-magnitude diagram (CMD) for all the ob- 
jects is presented in Fig. ^. As before, the objects retained 
as good candidates after visual inspection of the images 
are plotted as filled symbols and the rest as open symbols. 
The objects from Field 1 are plotted as diamonds, the ones 
from Field 2 as triangles. Vertical lines indicate the color 
range within which the MW and majority of M31 globular 
clusters are located. Note that several open symbols are 
found in the same area of the CMD as globular clusters. 
Spectroscopic observations are necessary to assess the real 
nature of these objects. 



Fig. 6. Color-Magnitude Diagram for globular clusters in 
NGC 5128 (filled circles) compared to the MW (stars) and 
M31 globular clusters (open squares). 

The object at (U-V)o = -0.99 and Vq = 19.96 
(fl.GC-8) is located very close to the area where a large 
number of bright young blue stars were found (Rejkuba 
et al. 2001a). It might be a young object, precursor of a 



globular cluster, similar to the ones found in other merg- 
ing galaxies (e.g. Whitmore et al. 1999), as well as to some 
clusters found in the central regions of NGC 5128 (Holland 
et al. |1999| ). 

There is a number of objects brighter than V < 17.5 
that were all identified in Field 1. Probably these are 
bright stars blended with some faint companions. They 
all appear saturated in the V-band VLT images and were 
found in the lower right corner (south-west) of Field 1 on 
U-band images. Their FWHM, that was larger than stel- 
lar, is probably an artifact due to the distortion of the 
PSF in that part of the chip. 

Figure I shows the CMD for aU NGC 5128 globular 
clusters (filled circles) that satisfy all the selection cri- 
teria. This CMD is compared with the ones of the MW 
(stars) and M31 (squares) globular clusters. Globular clus- 
ters populate a similar region of the CMD in the three 
galaxies. Data for the MW clusters were taken from the 
web page of W. Harris ( |1996D. M 31 clusters are from the 
catalogue of Barmby et al. ( 2000] ). 1 did not take into ac- 
count globular clusters with colors bluer than B— V= 0.55, 



as they might be younger objects (Barmby et al. 2000) 



The magnitudes were corrected for the mean reddening 
towards M31 (E(B-V) = 0.08). For the globular cluster 
candidates in NGC 5128 a mean for egrou nd reddening of 
E(B-V) = 0.1 (Burstein & Heiles |l982t Schlegel et al. 
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1998) was used. The reddening vector plotted in the up- 
per right corner corresponds to E(B— V) — 0.2. 

By constraining the colors to be within the range of 
the reddest MW and bluest M31 globular cluster colors 
(cf. the third selection criterion) I exclude distant compact 
galaxies (Barrientos ^000 ). The total number of bona fide 
globular cluster candidates is 71, with 23 in Field 1 and 
48 in Field 2. Only those are taken into account in the 
further analysis. 

Ks band magnitudes were measured on SOFI images 
for 23 of the bona fide globular cluster candidates and 
10 of the extended objects, some of which might also be 
globular clusters. The U— V color is strongly dependent 



on both metallicity and age (Buzzoni 1995; Kurth et al 



1999; Maraston 1998). Adding the K-band magnitude in- 



formation it is possible to disentangle the young from 
the old objects. The color-color plot is shown in Fig. 0. 
Filled dots are used to denote globular cluster candidates 
that satisfied all of the classification criteria. The lines are 
isochrones of evolutionary synthesis models of simple stel- 
lar populations (SSP) calculated for a Salpeter IMF and 
a range of metallicities from Z=0.0001 to 0.05 (see cap- 
tion for the exact metallicities plotted for each model). 
The four different models displa yed on panels (a), (b), (c) 
and (d) a re from Kur th et al. ( |l999|) , Br uzual & Char- 
lot (|200lD , Maraston ( ^00 1| ) and Worthey ( |l994| ), respec- 
tively The ages are indicated with numbers denoting Gyr 
next to the lines. Unfortunately the models do not allow 
one to discriminate between ages in the range of 5-20 Gyr. 
However, it is comforting to see that almost all the objects 
are concentrated towards the isochrones in excess of 1 Gyr. 
This is additional evidence that they are indeed globular 
clusters. A spread of metallicities of more than 1.5 dex is 
also evident. 

In Fig. most of the clusters lie slightly to the right 
and below the model lines. There are two possible expla- 
nations. First, a difference with the SSP models, particu- 
larly significant in U-V and B-V colors, has already been 
pointed out by Barmby & Huchra ( 200C ) . They suggested 
U fluxes from the stellar libraries, systematic problems 
with modeling HB color and observational errors as possi- 
ble causes of this offset. Second, in the data presented here 
additional offset may arise from the difference between the 
Bessel and Johnson U-band transmission curves (Bessel 



1995) 



4. Luminosity Function 

The V-band GCLF for NGC 5128 is compared to the ones 
of the MW and M31 in Fig. |. The actual shape of the 
luminosity function depends on the size of the bins and on 
the bin-centers. Statistically it is possible to compare them 
independently of binning using the Kolmogorov-Smirnov 
(KS) test. In Fig. || the cumulative distribution functions 
for the three GCLFs are shown. The thick line is used for 
NGC 5128, thin hne for the MW and dotted line for M31. 



Absolute magnitudes for NGC 5128 globular clusters were 
calculated using the distance modulus value of (m— M)o = 
27.8 (Soria et al. |1996D and reddening value of E(B-V) = 
0.1. 

The KS statistic measures the maximum value of the 
absolute difference between two cumulative distribution 
functions from which a probability that the two distri- 
butions are drawn from a common distribution can be 
calculated. The probability that the GCLF of NGC 5128 
matches the one of the MW is 89%, comparable to the 
probability that M31 GCLF matches the MW one (85%). 
The difference is slightly bigger in the U-band distribu- 
tions (Fig. fright panel), but it is also probable that the 
three luminosity functions are essentially similar for the 
three galaxies. 1 conclude that the NGC 5128 GCLF is as 
similar (or as different) to the MW GCLF as the M31 one. 

Usually the GCLF is fitted with a gaussian or ^5 dis- 
tribution function in order to determine its peak magni- 
tude. It should be noted, however, that the reasonably 
complete luminosity functions are known only for the two 
spiral galaxies, the MW and M31. Both of these GCLFs 
appear to depart from the gaussian distribution (Ashman 
et al. 1995), while there is no complete GCLF for any ellip- 
tical galaxy. The dependence of the turn-over magnitude 
on the detection of the faint end of GCLF is discussed in 
Ashman et al. (1995). The determination of the turn-over 
magnitude for NGC 5128 is left for a subsequent paper 
where the complete analysis of WFl data is going to be 
presented. 

The detection of the faint end of the GCLF in an ellip- 
tical galaxy is also important from the dynamical point of 
view. The faint globular clusters are typically also the less 
massive ones and their presence in the deep potential of 
a giant elliptical puts strong constraints on the potential 
itself as well as on the effects of tidal forces in the halo. 
Unfortunately the faintest globular clusters, correspond- 
ing to the MW Pal 1, AM 4, Terzan 1, Pal 13 are probably 



not detectable on FORSl images (see Sect. 2.4). 

The processes that destroy globular clusters can be 
divided into internal, like two-body relaxation and mass 
loss during stellar evolution, and external. The external 
processes depend greatly on environment and therefore 
on the position of the clusters in the galaxy. While tidal 
shocks preferentially destroy less dense clusters, dynamical 
friction brings more massive clusters very close to galactic 
centre. Both of these external processes are much stronger 
in the inner part of the galaxy. According to theoretical 



simulations (Vesperini 2000, Ostriker & Gnedin 1997) in 
a galaxy hkeNGC 5128, with the total mass of ^ 4 x 10" 
Mq (Israel |1998D and Re = 5.24 kpc (Dufour et al. |1979D, 



a strong disruption and mass loss are efficient only within 
a galactocentric distance of 1 — 2 Re, while at larger radii 
most of the mass loss is due to effects of stellar evolution. 
Field 1 is located at 2 - 3.5 Re and Field 2 at 1.2 - 2.3 
Re . While the dynamical effects may still be important for 
Field 2 clusters, they are probably completely negligible 
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Fig. 7. Color-color diagram for globular clusters in NGC 5128 (filled circles). The lines are the isochrones of SSP models. 
The number next to each isochronc indicate the age in Gyr. Each isochron e is pl otted for a range of metallicitics and 
for a Salpeter IMF. The four different mode ls are from: (a) Kurth et al. (|l999| ) for Z=0.0001, 0.0004, 0.008, 0.004, 
0.02 a nd 0.05; (b) Bruzual & Chariot ( 2001 ) for the same range of age and Z as the previous model; (c) Ma raston 
(|200l|) for [Fe/H]= -2.25 (only 14 Gyr isochrone), -1.35, -0.55, -0.33, 0.00 and 0.35 dex; and (d) Worthey (|l994|) 
for Z=0.0002, 0.0004, 0.008, 0.004, 0.02 and 0.05. 



for clusters in the outer field. It is not possible to quantify 
these effects by simple comparison of the two subsamples, 
because of the small number statistics. 



5. Color and Metallicity Distributions 

The (U— V)o color histogram for 71 globular clusters in 
NGC 5128 indicates a bimodal distribution (Fig. |lO|). Note 
that the (U— V)o color distribution of the MW globular 
clusters does not show bimodality while the one of M31 
does. The (U— V)o color distribution of the MW globular 
clusters is in fact biased due to lack of U-magnitude and 
higher reddening values (only clusters with E(B— V) < 0.5 
are plotted) of bulge (more metal-rich and redder) cluster. 
The blue peaks of the MW, M31 and NGC 5128 globular 
cluster color distributions coincide. The blue cutoff in the 
three galaxies is quite sharp and is driven principally by 



the colors of blue horizontal branch stars. The red side of 
the color distribution is much broader. In M31 it can be 
partially driven by differential reddening inside the host 
spiral galaxy, because the colors for M31 clusters presented 
here were corrected only for the mean foreground redden- 
ing of E(B-V)=0.08. Examination of the stellar (U-V) 
vs. (V— K) diagram in the two fields in NGC 5128 indi- 
cates that the amount of reddening does not vary much 
in the halo of this giant elliptical (Rejkuba et al. 2001b). 
Harris et al. (1992) and Holland et al. (1999) also found 
more red (metal-rich) clusters in their NGC 5128 data. 

The color distribution of the MW globular clusters in 
Fig. |l^, being dominated by halo clusters, is representative 
of the GCS of our Galaxy as if observed from outside. The 
NGC 5128 clusters in Harris et al. ( |1992D and in this work 
belong as well to the halo population. The differences in 
the color (metallicity) distribution of clusters and stars in 
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Fig. 8. From top to bottom: V-band GCLF for M31, 
NGC 5128 and MW. 




Fig. 9. Kolmogorov-Smirnov statistic: The cumulative 
distribution function for NGC 5128 (thick line) is con- 
structed from the V-band GCLF (left panel) and U-band 
GCLF (right panel) and compared with the ones of the 
MW (thin line) and M31 (dotted line). 



the halo of a spiral, like MW and a gE hke NGC 5128, 
are reflected in figure 10 with a much larger metallicity 
spread in gE galaxy. 

The line overplotted on the color distribution of 
NGC 5128 is the best nonparametric kernel estimate (Sil- 
1986| ) of the color distribution used to measure 



verman 



accurately the blue and the red peaks of the distribu- 
tion. The blue and red peaks are measured at (U— V)o 
0.73 ± 0.01 and 1.28 ± 0.01, respectively. 



From the color-color diagram (Fig. |^) it is obvious that 
the globular clusters in NGC 5128 are older than 1 Gyr. 
Unfortunately, the models do not allow one to better con- 
strain the ages (Fig. |^; see also Barmby & Huchra 2000). 
Assuming that the globular clusters in NGC 5128 have 
similar ages as the MW ones, one can use the Galactic 
globular clusters to calibrate the metallicity scale. It is 
known that for the very low and very high metallicities the 
color-metallicity relation cannot be expressed by a simple 
linear function, but extrapolation from lower metallicities 
towards the higher ones by using the higher polynomial 
fit might lead to unphysical results. Thus I prefered to 
make a linear fit with gaussian errors in a^Fe/H) — 0.1 
in order to calibrate the metallicity as a function of in- 
trinsic ( U— V )o color. Only MW globular clusters from W. 
Harris ( |l996[ ) web list with E(B- y)< 0.5 (equivalent to 
E(U-V)< 0.82; Rieke & Lebofsky |l985| ) were used in the 
fit. The follow ing co lor-metallicity relation is based on the 
Zinn & West (1984) scale and is valid only for the range of 
metallicities observed in the MW globular clusters (Fig. 
0): 

[Fe/H] = -2.93(±0.07) + 1.79(±0.07) x {U - V)o (4) 



In the Fig. ^ hnear (U-V)o vs. [Fe/H] color- 
metallicity relation (Eq. ^; full line) is plotted along with 
the evolutionary synthesis models of simple stellar pop- 
ulations (short dashed lines) from Kurth et al. (1999). 



Using the empirical linear calibration and the SSP mod- 
els I calculate the metallicity of the two peaks in the 
color-distribution. The blue and the red peak correspond 
to [Fe/H]= -1.7 dex and [Fe/H]= -0.6 dex using the 
linear relation (Eq. ^) and to [Fe/H]= —1.7 dex and 
[Fe/H]= —0.5 dex using the models. 
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Although it is not obvious from the (U— V)o color 
distribution, but rather from the spectroscopically mea- 
sured [Fe/H] distribution, the metallicity distribution of 
MW globular clusters is bimodal with the metal-poor and 
metal-rich peaks at -1.6 and -0.6 dex respectively (Har- 
ris 200C| ) . These values are practically the same as the ones 
derived here for globular clusters in NGC 5128 from the 
(U— V)o colors. On the other hand, on the basis of Wash- 
ington photometry of 62 globular clusters in NGC 5128, 
Harris et al. ( 1992 ) conclude that "in both the range and 
distribution of abundance, the NGC 5128 GCS is clearly 
different from the Galactic GGS and quite similar to that 
of the large elliptical NGC 1399". They measured a mean 
metallicity of ([Fe/H](^_j,^) = —0.8 ± 0.2 dex, substan- 
tially more metal-rich than the mean value for the MW 
GCS. They, however, mention that the values exceeding 
[Fc/H]'^ —0.25 are extrapolated, a fact that is also true 
in the metallicity calibration from (U— V)o color distribu- 
tion. The extrapolation leads to large errors and a possible 
overestimation of abundances of the metal-rich clusters. 
Because of this I do not report here the metallicities of 
individual clusters. 

Note also that I restricted the (U— V)o color to be 
within the range observed in MW and bulk of the M31 
globular clusters. In this way the sample was selected to be 
as clean as possible from the very compact and red galax- 
ies. It is not surprising that the clusters with supra-solar 
abundances are not present, reducing the mean abundance 
of the sample. However, relaxing the red color cut and 
adding the 6 "good" globular cluster candidates redder 
than (U-V)o = 2.5 (filled symbols in Fig. |) the mean of 
the color distribution shifts from ((U - V)o) = 1.05 ±0.45 
to ((U — V)o) — 1.22±0.75, corresponding to mean metal- 
licities of ([Fe/H]) = -1.1 ± 0.1 dex and ([Fe/H]) = 
—0.7 ±0.2 dex, respectively. The latter value is very simi- 
lar to the mean metallici ty der ived from Washington pho- 
tometry by Harris et al. ( 1992 ). In order to assess the real 
nature of the objects redder than (U— V)o = 2.5 spectro- 
scopic observations are necessary. 

6. Conclusions 

The main results and their implications presented in this 
work are the following: 

1. On the basis of high resolution ground-based images, 
taken with VLT and FORSl, I identified 71 bona fide 
globular clusters in the halo of the nearest giant ellip- 
tical galaxy NGC 5128. Only 5 of these clusters were 
previously known. 

2. In a UV color-magnitude diagram for the complete 
sample of 71 clusters and a UVK color-color diagram 
for a subsample of 23 clusters, the objects span a simi- 
lar magnitude and color range as the globular clusters 
in the MW and M31. Note, however, that there is no 
cluster corresponding to uj Gen, the brightest MW clus- 
ter, in these two fields. The position of 23 objects in the 
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Fig. 11. The (U— V)o vs. [Fe/H] color-metallicity linear 
relation (full line) derived from the MW globular clusters 
with E(B-V)< 0.5 (stars). The M31 clusters with spec- 
troscopic metallicities and accurate (U-V) colors plotted 
for comparison (squares) were not used for the fit. Typical 
errors in measured metallicities are plotted on in the box. 
Short dashed lines are the SSP models from Kurth et al. 
( 1999| ) for ages of 5, 10 and 15 Gyrs, from blue to left, re- 
spectively. Each model spans metallicities from Z=0.0001 
to 0.05. 



(U— V)o vs. (V— K)o color-color diagram is consistent 
with their classification as old globular clusters. 
The GCLFs spanning -10.1 < My < -4.9 and 
—9.3 < AIjj < —3.3 have been constructed. These 
are the deepest GCLFs of an elliptical galaxy made 
so far. Kolmogorov-Smirnov statistics show that the 
difference between the GCLFs of NGC 5128 and MW 
is not larger than the difference between the GCLFs of 
M31 and MW. Similarity of the GCLFs of an elliptical 
with respect to the spiral galaxy had never been tested 
at the faint end before. 

The presence of faint globular clusters in the halo of 
NGC 5128 puts constraints on the effectiveness of the 
tidal forces in the deep elliptical galaxy potential. The 
dynamical effects may be important for clusters that 
are found within ~ 2 Re from the galactic center. Un- 
fortunately, the selection of globular clusters from VLT 
images is not sensitive to the most compact and the 
faintest clusters, similar to ones like Pal 1, Pal 13, 
AM 4 and Terzan 1 in our Galaxy, while the less dense 
ones start to be confused with background galaxies at 
faint magnitudes. 

The (U— V)o color histogram of 71 clusters indicates 
a bimodal distribution, supporting the Zepf & Ash- 



man (1993) suggestion. Assuming that the clusters in 



NGC 5128 and in the MW span a similar age inter- 
val, and adopting the linear fit between the (U— V)o 
color and metallicity, I derived the [Fe/H] of the red 
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and blue peaks of the bimodal distribution to be —1.7 
dex and —0.6 dex, respectively. Using the SSP mod- 
els from Kurth et al. (1999), instead of the linear 
fit, the results do not change significantly, giving val- 
ues of —1.7 dex and —0.5 dex for the metal-poor and 
metal-rich peaks, respectively. This is different from 
the Harris et al. ( |1992D and Zepf & Ashman ( [1993D 
result, but is partially due to the adopted color cut 
of (U— V)o < 2.5 in selecting globular clusters and 
to small number statistics, since both samples have 
< 10% of the total cluster population. Relaxing the 
red color cut the mean color of the distribution cor- 
responds to mean metallici ty of [Fe/H] = —0.7 dex, 
very similar to Harris et al. ( 1992| ) result. The true na- 
ture of the redder objects has to be assessed through 
spectroscopy. 
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